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Abstract 

We report recent results on the role of instabilities in the isotropization and thermalization of a longitudinally 
expanding system of quantum fields. 
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1. Introduction 

The evolution of the quark-gluon matter produced in ultrarelativistic heavy ion collisions is apparently well 
described by hydrodynamics with a very small viscosity. Given its derivation from an expansion of the energy- 
momentum tensor in gradients of the velocity field, one expects hydrodynamics to be applicable only when the trans- 
verse and longitudinal pressures are sufficiently close, and if the system is close enough to local thermal equilibrium. 

An outstanding problem in the application of Quantum ChromoDynamics to the description of the early stages 
of these collisions is to justify that these conditions are indeed met at a short enough time to apply hydrodynamics. 
The main QCD tool for performing these studies is the Color Glass Condensate (CGC) framework [1 1, an effective 
theory designed to describe the highly occupied gluon states encountered in the wave functions of high energy nuclei, 
and their interactions in a collision. In the CGC, the fast partons (k + > A + for a nucleus moving in the +z direction) 
are described as static (thanks to time dilation of all the internal dynamics in a high energy nucleus) color sources 
jf ~ fif+p on the light-cone, with a probability W[p]. The gluons below the cutoff (k + < A + ) are described as ordinary 
gauge fields, coupled eikonaly to the fast sources. When this effective description is applied to a collision, the current 
y is the sum of two terms, corresponding respectively to each of the colliding nuclei. 

2. Initial state factorization 

An important issue in applications of the CGC to the calculation of observable quantities is the dependence on 
the cutoffs that are introduced to separate the fast and slow partons. This cutoff is not a physical quantity and should 
cancel in observables. However, since it is the upper bound in loop integrals, higher order corrections usually produce 
logarithms of the cutoff. It is well known that in electron-nucleus collisions, these logarithms can be absorbed by 
letting the distribution W[p] depend on the scale A + according to the JIMWLK equation [2j- For the CGC to provide 
a consistent description of nucleus-nucleus collisions, one needs to prove that the same distributions W[p]'s are 
sufficient to absorb all the logarithms that arise in these collisions, and that this works for a sufficiently large class of 
observables. 

The first step in this proof is to show that at Leading Order (LO) in the strong coupling a s , all the inclusive 
observables (i.e. those that do not restrict the content of the final state, such as particle spectra, or the expectation 
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value of some local operators in the final state) can be expressed in terms of the retarded solution of the classical 
Yang-Mills equations that vanishes in the remote past (x° — > -oo) [3 1. For instance, in terms of this classical field JV, 
the gluon spectra at LO read 
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The next step is to compute the Next to Leading Order (NLO) correction to these observables. In 0, we have 
shown that for an inclusive observable O, the LO and NLO contributions are formally related by 



a, 



(2) 



In this formula, T is the operator that generates shifts of the classical field at some initial time (and the integrals over 
u, v are integrals over the spatial coordinates at this time). The quantities a are small perturbations to the classical 
field, evaluated at the same time. The formula is true for any initial time, but some choices simplify the problem by 
making a^, a calculable analytically. Then, the main result in proving the factorization of the logarithms of the cutoffs 
is 

\f f[a k T]„ K T] v + j [or T] B = log (A + ) <Hi + log (A") <H 2 + terms w/o logs , (3) 

where "Hi^ are the JIMWLK Hamiltonians of the two nuclei. This formula shows that the logarithms of the two cutoffs 
do not mix, since they are multiplied by operators that act on a single nucleus. From this formula, it is easy to show 
that one can absorb all the leading logarithms in two distributions W^2 that evolve with their respective JIMWLK 
equation, 

(O)^ L „ g = j [DP* DP*] Wi [/>,] W 2 [p 2 ] LO [p uP2 ] . (4) 

It is the universality of the distributions W (i.e. the fact that the same distributions enter in different types of reactions 
involving the same nuclei) that gives its predictive power to the CGC. 



3. Thermalization and isotropization 



3.1. Resummation of the unstable modes 

At LO, the classical color field produced in a nucleus -nucleus collisions corresponds to longitudinal chromo- 
electric and chromo-magnetic fields (figure [T|- left) that form flux tubes whose transverse size is of the order of the 
inverse saturation momentum. The energy momentum tensor of this field configuration has a negative longitudinal 
pressure, and is thus quite far from what can be reasonably described by hydrodynamics. 

Moreover, it has been known for some time that the classical solutions of Yang-Mills equations have instabilities: 
small rapidity dependent perturbations to the boost invariant classical field encountered at LO grow exponentially in 
time, and eventually become as large as the classical field itself [5]. Such instabilities occur in loop corrections and 
lead to secular divergences in observables. These unstable contributions can be tracked by an improvement of the 
power counting, where one assigns an exponential of time to each operator T. Then one sees that the terms with the 
strongest divergences (one exponential factor per power of the coupling g - see figure [T] right, for examples of 2-loop 
graphs that differ in the strength of their secular divergences) are obtained by exponentiating the operator in eq. pi. 
This exponentiation performs a resummation of a selected subset of all the higher loop corrections. When applied to 
the energy-momentum tensor, this resummation reads 



- exp[iJ^ M „KT] v + J[aT]„ 
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Figure 1: Left: color flux tubes. Right: examples of 2-loop corrections. 

The second line is an identity, based on the fact that the operator T is the generator of the shifts of the initial classical 
field ytinit. It shows that this resummation amounts to a Gaussian average over fluctuating initial conditions, centered 
on the classical one. 

3.2. Non expanding system 

In order to test the effect of this resummation in a simpler setting, we have studied a scalar field theory with 
a 4 interaction. This theory has also unstable classical solutions, due to parametric resonance. Moreover, in four 
space-time dimensions, this theory is scale invariant at the classical level, as is Yang-Mills theory. 

We have first considered this theory in the case of a fixed volume @. First of all, one can check explicitly the 
presence of secular divergences in fixed loop order computations of the energy-momentum tensor (see figure [2] left). 
On the right plot of this figure, we also see that the behavior is quite different after the resummation: not only has the 
secular divergence been cured, but one also sees that the pressure relaxes to the equilibrium value P = e/3. 




Figure 2: Left: energy density and pressure at LO and NLO. Right: resummed energy density and pressure. 

In this set-up, one can also compute the occupation number at various times in the evolution (figure [3] left). In 
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Figure 3: Left: evolution of the particle distribution. Right: time evolution of the occupation number at k = 0. 
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this computation, the initial classical field was set so that the occupation number is large in a single momentum 
mode, and zero elsewhere. Rapidly, a population of particles builds up in the neighboring modes. At late times, the 
distribution has evolved into a classical equilibrium. Interestingly, the best fit requires a positive chemical potential. 
This indicates that chemical equilibration has not yet happened at these times, suggesting that the number changing 
processes are considerably slower than the elastic ones. Moreover, this chemical potential is almost equal to the 
mass of the quasi-particles, and sees a marked excess of particles in the zero mode, both indicative of Bose-Einstein 
condensation (BEC). By following the occupation of the zero mode over longer times (figure [3] right), one sees that 
indeed this BEC is only transient and disappears eventually when inelastic processes equilibrate the particle number. 
Whether a similar phenomenon can occur for gauge fields is still an open issue |7). 

3.3. Longitudinal expansion 




Figure 4: Left: relaxation of the trace of T^ v . Right: isotropization of the pressures. 

The previous scalar model can be studied also for a longitudinally expanding system, where the classical back- 
ground field is boost invariant [8 |. In the left panel of figure |4] one first sees that the trace of the pressure tensor 
relaxes towards the energy density. Moreover, the time dependence of the energy density changes from t _1 to r~ 4 / 3 , 
indicating a qualitative change from a very small longitudinal pressure to a longitudinal pressure equal to e/3. This is 
confirmed in the right panel of figure |4] after a very rapid initial decrease, due to the redshifting of the longitudinal 
momenta, the longitudinal pressure starts growing again due to the instabilities, to eventually become very close to 
the transverse pressure. The spectrum of fluctuations appropriate for gauge fields over a boost invariant background 
has been worked out in [9|, paving the way for the same simulation in Yang -Mills theory. 
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